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ABSTRACT 

We present a deep broadband optical imaging study of a complete sample of luminous elliptical 
galaxies {Mb < —20) at distances 15 Mpc - 50 Mpc, selected from the Tully catalog of nearby galaxies. 
The images are flat to ~0.35% across the 20' field and reach a V band depth of 27.7 mag arcsec"^. 
We derive an objective tidal interaction parameter for all galaxies and find that 73% of them show 
tidal disturbance signatures in their stellar bodies. This is the first time that such an analysis is done 
on a statistically complete sample and it confirms that tidal features in ellipticals are common even 
in the local Universe. From the dynamical time of the sample galaxies at the innermost radius where 
tidal features are detected we estimate the mass assembly rate of nearby ellipticals to be dM/M ~ 0.2 
per Gyr with large uncertainty. We explore the relation between gravitational interaction signatures 
and the galaxy environment and find that galaxies in clusters are less disturbed than group and field 
galaxies. We also study how these interactions affect the broadband colors of ellipticals and find a 
moderate correlation, suggesting that the mergers are not accompanied by significant star-formation. 
Lastly, we find no correlation between AGN activity, as measured by 6cm radio emission, and large 
scale tidal distortions. This implies that gravitational interactions are not the only, and perhaps not 
the most important, trigger of nuclear activity. In summary, we find that elliptical galaxies in groups 
and low density environments continue to grow at the present day through mostly "dry" mergers 
involving little star formation. 

Subject headings: galaxies: interactions — galaxies: evolution — galaxies: elliptical — galaxies: 
structure 



1. INTRODUCTION 

Giant elliptical galaxies dominate the high end of the 
mass spectrum in the nearby Universe and are mostly 
composed of old stellar populations, in contradiction to 
naive predictions from hierarchical growth models. In 
a ACDM Universe the most massive dark matter ha- 
los form later than less massive ones and the elliptical 
galaxies that reside in them are therefore expected to 
have formed recently (e.g., Kauffmann et al. 1993). This, 
however, does not agree with photometric studies of el- 
lipticals that produce red colors in all of these galaxies, 
implying their stars are mostly old (e.g. Visvanathan & 
Sandage 1977; Bower et al. 1992; Chang et al. 2006). 
The lack of recent star formation is further supported by 
results from HI surveys of early-type galaxies that find 
extremely low cold gas content in these systems when 
compared to spiral galaxies (e.g. Sanders 1980; Knapp 
et al. 1985; Sadler 2001). 

A possible solution may be that galaxies grow through 
"dry" mergers which do not lead to star formation (Kor- 
mendy 1984; White & Frenk 1991; Kauffmann 1996; van 
Dokkum 2005; Naab et al. 2006; Boylan-Kolchin et al. 
2006). The little gas that is accreted into the elliptical 
by the interaction may be heated and kept from cooling 
by nuclear feedback (Croton et al. 2006) or the gravita- 
tional interactions themselves (Dekel & Birnboim 2008; 
Kenney et al. 2008). 

One of the most direct ways to study the assembly 
rate of elliptical galaxies is to measure the incidence of 
tidal features. Schweizer & Seitzer (1992) studied a sam- 
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pie of 69 early-type galaxies with a recession velocity of 
less than 4000 kms~^. The galaxies were selected from 
a parent sample of 145 objects to exclude most cluster 
members and to include a few galaxies that had been 
previously known as interacting systems. Optical colors 
were measured from the sample and a "fine structure" 
parameter was assigned to each galaxy by visually char- 
acterizing the amount and type of their morphological 
disturbances. In their paper, Schweizer & Seitzer found 
that the broadband colors of early-type galaxies are cor- 
related with their respective fine structure parameter and 
therefore also with gravitational interactions. 

A different approach to quantifying the morphologi- 
cal disturbances of the stellar bodies of ellipticals was 
performed by van Dokkum (2005) [hereafter vD05], who 
studied a sample of 126 early- type galaxies at a median 
redshift of 0.11. vD05 divided deep broad-band optical 
images of the galaxies by their respective fitted smooth 
stellar models and measured the residuals of the quo- 
tients. This quantitative approach facilitates comparison 
to simulations (e.g. Feldmann et al. 2008) and other data 
(e.g. Mcintosh et al. 2008). 

Here we perform, for the first time, a quantitative anal- 
ysis of morphological disturbances in a complete sample 
of elliptical galaxies at z?^0. The field of view of current 
generation imagers, along with excellent flat fielding, al- 
low us to study the stellar bodies of nearby ellipticals and 
to constrain the rate of gravitational interaction in the 
nearby Universe. The results of this study can be used in 
numerical simulations and theoretical models of galactic 
evolution to constrain the formation and mass growth of 
these objects at recent cosmic times. Furthermore, most 
of our sample is composed of well-known bright galaxies 
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Fig. 1. — The luminosity and environment distributions of the 
55 sample galaxies. 



that have been studied extensively in the past. 

2. DATA 

Observations of a volume limited sample of luminous 
elliptical galaxies were obtained with the Y4KCam im- 
ager mounted on the SMARTS Im telescope at CTIO in 
four epochs between 2006 and 2008. The instrument is a 
4Kx4K CCD camera optimized for wide-field broad-band 
imaging, providing a nearly undistorted 20'x20' field of 
view. 

2.1. Sample Selection 

Candidates for inclusion in this study were initially se- 
lected from the Nearby Galaxies Catalog (Tully 1988), 
consisting of all elliptical galaxies at a declination be- 
tween -85 and +10. A distance cut was applied to the 
initial sample to exclude galaxies that are farther than 
50 Mpc or closer than 15 Mpc. The distance thresh- 
old ensures that the outer parts of all galaxies fit in the 
instrument's field-of-view and that a sufficient signal-to- 
noise ratio is achieved for all targets. A luminosity cut of 
Mb <-20.15 was used, with the magnitudes taken from 
Tully (1988) and converted to our cosmology. This limit 
corresponds to Mb <-20.0 in the Tully (1988) atlas, as 
he used Hq = 75 km s~^ Mpc~^, and to L>L^ (Blanton 
et al. 2003). The last selection criterion that we used 
rejected galaxies with Galactic latitude of less than 17 
due to difficulties of constraining a good model fit in a 
crowded stellar field. 

One galaxy (NGC 5796) was excluded due to poor ob- 
serving conditions that resulted in noise levels over 5 
times worse than the sample mean. The final catalog 
therefore consists of 55 giant ellipticals (figure 1) and it 
includes members of four nearby clusters (Virgo, Fornax, 
Centaurus and Antlia). Galaxy environments were de- 
termined from the literature using NASA's Astrophysics 
Data System Bibliographic Services to find references for 
previous studies of the objects. The methods used to de- 
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Fig. 2. — The observing pattern of a galaxy at (0,0) overlaid 
on an SDSS image of NGC 4636. Filled circles represent telescope 
pointings for both target and sky frames, which are accompanied 
by a number denoting their place in the sequence. The central 
square shows the angular size of the field of view. 



rive these assumed environments are therefore inconsis- 
tent throughout the sample and they vary in accuracy. 

2.2. Observations 

Broadband optical observations of the sample galaxies 
had been acquired by multiple authors, showing some of 
the tidal features that we discuss in this work. However, 
the standard techniques used to obtain and reduce these 
previous data were not aimed at revealing the faint gravi- 
tational interaction signatures and they were washed out 
by residual background level variations across the field. 
In order to overcome this we supplemented each set of 
galaxy observations with a sequence of dark sky expo- 
sures of equivalent depth. 

The galaxies were observed in a sequence of 33 point- 
ings, including both object and dark sky frames. We 
used the V band as it provides the highest signal-to-noise 
ratio in a given exposure time. The observing pattern, 
shown in figure 2, included seventeen 300 sec exposures 
of the target galaxy and sixteen 300 sec exposures of 
background sky, acquired in alternating order. All ob- 
ject frames were shifted from each other by 1-3' in order 
to correct for cosmetic defects in the CCD. Figure 2 also 
shows the order in which sky frames were observed, col- 
lecting data from different sides of the galaxy in a semi- 
random pattern. All galaxies were observed with total 
exposure times between 4200 and 7200 seconds. 

In order to improve the signal-to-noise ratio the data 
were binned by a factor 2x2 at the telescope, producing 
a pixel size of 0.578". To further increase the sensitivity 
to tidal features we also binned in software, resulting in 
a final pixel size of 1.156". The typical stellar FWHM of 
the images used for analysis is ~1.7". 

2.3. Reduction 

Initial reduction of the data followed standard tech- 
niques and consisted of zero level subtraction and first 
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Fig. 3. — A demonstration of the importance of a high degree of field flatness for studying faint interaction signatures. The same data 
were reduced with and without the use of a dark-sky flat field frame (left and right panels, respectively). Note that in the left panel 
the tidal tail extending away from NGC 5576 cannot be distinguished from fiattening distortions in the field. The three boxes in the 
upper-right corner have surface brightness values similar to the faintest detectable features. Both images were further smoothed and binned 
for presentation purposes. 



order field variation corrections using dome fiat frames. 

In order to flatten the field to a higher degree dark- 
sky fiat-field frames were prepared for each of the sam- 
ple galaxies. All sources were masked in the individual 
dark-sky exposures which were then averaged and ap- 
plied to the object frames. Lastly, the reduced object 
exposures were aligned and combined to create the final 
data products. 

Apparent magnitudes were calibrated using aperture 
photometry of Prugniel & Heraudeau (1998) and were 
corrected for Galactic reddening using infrared dust 
maps from Schlegel et al. (1998). We assume distance 
measurements from the Tully catalog (corrected to our 
cosmology) to convert the luminosity profiles to physical 
units. 

The data were reduced using the NOAO/IRAF 
software and are available for download at 
http:/ /www. astro.yale.edu/obey 

2.4. Depth and Flatness of the Field 

Our sky fiat-fielding procedure (discussed in subsection 
2.2) resulted in a dramatic improvement to the flatness 
of the field, as can be seen in figure 3. In this figure 
the ~85 kpc long tidal tail that extends away from the 
companion of NGC 5576 is barely detectable in the image 
on the left which uses standard fiat-fielding but is clearly 
visible in the dark-sky fiat-fielded frame. 

We used two methods to determine the depths of the 
sky fiat field images. First, we reduced and stacked a 
dark sky image in the same way that was used to pro- 
duce the target frames and smooth it to the scale of a 
typical tidal feature (~20"). Although the resultant im- 
age was inherently fiat, it preserves the photon noise level 
as limited by the telescope, instrument and the observing 
program. Second, we aggressively masked a target frame 
for bright objects and subsequently smooth it using a 



median kernel. This frame is insensitive to any pixel-to- 
pixel variation but it reflects the large-scale variation due 
to residual fiattening issues or sky conditions at the time 
of observation. We then measured the standard devia- 
tion across both frames to obtain the detection threshold 
of faint tidal features. For a 20" x20" box (corresponding 
to a typical tidal feature scale) we derived a la photon- 
noise detection threshold of ~29 mag. From the flatness 
limited frame we measured a la detection threshold of 
~27.7 mag. This result implies that a similar program 
carried on a larger telescope will not necessarily yield a 
lower detection threshold as the data are dominated by 
residual flatness variations rather than photon noise. 

3. TIDAL FEATURES 
3.1. Morphological Disturbances 

A visual inspection of the 55 sample galaxies reveals 
an extraordinary assortment of stellar morphology dis- 
turbances in a large number of galaxies. Similar distur- 
bances have been reported by several authors and they 
can be divided into four main categories: 

Shells - observed in the stellar component of many 
nearby ellipticals, shells are usually considered to result 
from an accretion of a small companion by a massive 
galaxy. Numerical simulations of such interactions have 
successfully reproduced the observed properties of shells 
even without the inclusion of gas dynamics and star for- 
mation (e.g. Dupraz & Combes 1986). In our sample we 
identify at least 12 galaxies with shells, making them the 
most common interaction signature. 

Tidal Tails - these linear streams of stellar matter 
are less commonly observed than shells in elliptical sys- 
tems and are generally accepted as evidence for a dy- 
namically cold component in the accreted companion. 
Tidal tails can also be reproduced in simulations of dry 
mergers when one of the progenitors is rotating (Combes 
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Fig. 4. — Three examples of the model fitting and division method we employ to derive the tidal parameter of the sample galaxies. From 
top to bottom: the galaxy NGC 5638 and its companion show no signs of gravitational interaction; NGC 3923 has shells at multiple radii; 
NGC 2865 is a recent or ongoing merger showing multiple interaction signatures. From left to right: The total counts image; foreground 
and background object masked frame; elliptical model fit to the data; the division of the data by the model - this image was used to 
measure the tidal parameter T. 



et al. 1995). Tails that are produced in this manner ap- 
pear wide and short-lived in comparison to the long and 
narrow tails produced by interacting spirals. An inter- 
action between an elliptical galaxy and a gas-rich disk 
system can also produce such features (Feldmann et al. 
2008). In this work we identify both kinds of tails in a 
range of widths and lengths, some extending as far as 11' 
(~85kpc) in projected length. 

Broad Fans of Stellar Light - also well reproduced 
in simulations of dry mergers, extended stellar fans typ- 
ically have a low surface brightness which makes them 
hard to detect in shallow surveys. Studies of interacting 
ellipticals (vD05, Mcintosh et al. 2008) show multiple 
examples of this morphological signature and associate 
them with interactions of gas poor systems. This inter- 
action signature is visible in a number of systems in our 
sample, most of which do not have an obvious compan- 
ion and are therefore likely the result of a major merger 
in the past (see figure 1 of vD05). 

Highly Disturbed Galaxies - at least three galaxies 
(^^5%) in our sample of 55 exhibit signs of an ongo- 



ing merger that is violently disturbing their stellar 
components (NGC 2865, NGC 3640 and NGC 5018). 
These systems show signs of one or more of the above 
interaction signatures, including well defined linear 
filaments that extend radially from the galactic centers. 
Two of the three galaxies are recent merger remnants 
and in the third the interaction is still ongoing with a 
nearby companion. These complex systems probably 
result from mergers with multi-component objects such 
as SO or spiral galaxies. 

Morphological disturbances in the stellar bodies of the 
sample galaxies are common and are clearly apparent 
in at least two-thirds of the objects. The existence and 
frequency of these features imply that elliptical galaxies 
are still evolving in the nearby Universe through colli- 
sions and gravitational interactions that affect their stel- 
lar bodies. 

A brief description of the visible distortions of the sam- 
ple galaxies is given in the last column of table 1. 
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TABLE 1 
Catalog 



Object 




Dec.^ 


Mb ^ 


{B - vy 


Env.^ 


Tc ^ 


Tidal Feature Description 


NGC 584 


01:31:20.7 


-06:52:05 


-20.83 


0.92 


F 


0.076 


Non-spherical isophotes; Spiral companion 


NGC 596 


01:32:52.1 


-07:01:55 


-20.37 


0.90 


F 


0.110 


Pronounced outer shell 


NGC 720 


01:53:00.5 


-13:44:19 


-20.53 


0.96 


F 


0.079 


Boxy outer isophotes 


NGC 1199 


03:03:38.4 


-15:36:49 


-20.24 


0.97 


G 


0.067 


Intra-group diffuse emission 


NGC 1209 


03:06:03.0 


-15:36:41 


-20.51 


0.95 


G 


0.116 


X-like isophotal structure; NW linear feature 


NGC 1399 


03:38:29.1 


-35:27:03 


-20.49 


0.98 


C 


0.064 




NGC 1395 


03:38:29.8 


-23:01:40 


-20.44 


0.94 


F 


0.094 


NW perpendicular feature 


NGC 1407 


03:40:11.9 


-18:34:49 


-21.17 


0.93 


G 


0.083 




NGC 2865 


09:23:30.2 


-23:09:41 


-20.85 


0.78 


F 


0.193 


Multiple shells and tail-like features 


NGC 2974 


09:42:33.3 


-03:41:57 


-20.78 


0.95 


F 


0.110 


Multiple shells 


NGC 2986 


09:44:16.0 


-21:16:41 


-20.85 


0.99 


F 


0.045 


Possibly paired with disturbed spiral 


NGC 3078 


09:58:24.6 


-26:55:37 


-20.78 


0.97 


F 


0.103 




NGC 3258 


10:28:53.6 


-35:36:20 


-20.47 


0.92 


G 


0.123 


Possibly paired with low mass elliptical 


NGC 3268 


10:30:00.6 


-35:19:32 


-20.43 


0.96 


C 


0.087 


Possibly paired with disturbed spiral 


NGC 3557B 


11:09:32.1 


-37:20:59 


-20.28 


0.86 


N 


0.182 


Isophotal twisting; Asymmetric outer isophotes 


NGC 3557 


11:09:57.6 


-37:32:21 


-22.33 


0.87 


G 


0.111 


SW fan; Asymmetric outer isophotes 


NGC 3585 


11:13:17.1 


-26:45:18 


-21.08 


0.91 


F 


0.048 


Asymmetric outer isophotes 


NGC 3640 


11:21:06.8 


-h03:14:05 


-20.93 


0.92 


F 


0.142 


Highly disturbed stellar body 


NGC 3706 


11:29:44.4 


-36:23:29 


-21.29 


0.93 


U 


0.120 


Inner shell 


NGC 3904 


11:49:13.2 


-29:16:36 


-20.26 


0.94 


U 


0.108 


Outer shell 


NGC 3923 


11:51:01.8 


-28:48:22 


-21.38 


0.95 


U 


0.100 


Multiple outer shells 


NGC 3962 


11:54:40.1 


-13:58:30 


-20.88 


0.95 


F 


0.059 




NGC 4105 


12:06:40.8 


-29:45:37 


-20.51 


0.87 


F 


0.109 


Interacting with a tidally disturbed spiral 


NGC 4261 


12:19:23.2 


+05:49:31 


-21.56 


0.98 


C 


0.053 


NW tidal arm; Faint SE fan 


NGC 4365 


12:24:28.2 


+07:19:03 


-20.67 


0.97 


C 


0.070 


Faint SW fan 


IC 3370 


12:27:37.3 


-39:20:16 


-21.38 


0.89 


U 


0.192 


X-like isophotal structure; Broad N fan 


NGC 4472 


12:29:46.7 


+08:00:02 


-21.97 


0.97 


C 


0.000 




NGC 4636 


12:42:49.9 


+02:41:16 


-20.83 


0.93 


F 


0.066 




NGC 4645 


12:44:10.0 


-41:45:00 


-21.03 


0.95 


C 


0.000 




NGC 4697 


12:48:35.9 


-05:48:03 


-21.82 


0.92 


C 


0.091 


Non-spherical inner isophotes 


NGC 4696 


12:48:49.3 


-41:18:40 


-22.20 


0.94 


C 


0.075 


Faint outer shell 


NGC 4767 


12:53:52.9 


-39:42:52 


-21.28 


0.93 


C 


0.000 


Faint inner shell 


NGC 5011 


13:12:51.8 


-43:05:46 


-21.25 


0.89 


U 


0.077 




NGC 5018 


13:13:01.0 


-19:31:05 


-21.61 


0.85 


F 


0.184 


Highly disturbed; Multiple tidal tails and shells 


NGC 5044 


13:15:24.0 


-16:23:08 


-21.16 


0.98 


G 


0.041 




NGC 5061 


13:18:05.1 


-26:50:14 


-21.34 


0.85 


F 


0.104 


Pronounced tidal tail; Shell 


NGC 5077 


13:19:31.7 


-12:39:25 


-20.67 


0.98 


G 


0.061 




NGC 5576 


14:21:03.7 


+03:16:16 


-20.55 


0.88 


F 


0.122 


Interacting pair; tidal tail longer than 75 kpc 


NGC 5638 


14:29:40.4 


+03:14:00 


-21.27 


0.94 


F 


0.036 




NGC 5812 


15:00:55.7 


-07:27:26 


-20.73 


0.94 


F 


0.080 


Interacting with dwarf companion; Tidal tail 


NGC 5813 


15:01:11.2 


+01:42:07 


-20.92 


0.95 


G 


0.054 




NGC 5846 


15:06:29.3 


+01:36:20 


-21.31 


0.98 


G 


0.068 


Faint outer shells 


NGC 5898 


15:18:13.5 


-24:05:53 


-20.64 


0.92 


G 


0.114 


Three spiral arm-like tidal tails 


NGC 5903 


15:18:36.5 


-24:04:07 


-21.03 


0.89 


G 


0.075 


Possible dust extinction S of the center 


IC 4797 


18:56:29.7 


-54:18:21 


-20.90 


0.92 


G 


0.226 


Tidal tails; SE diffuse emission 


IC 4889 


19:45:15.1 


-54:20:39 


-20.70 


0.88 


U 


0.158 


Isophotal twisting; E fan 


NGC 6861 


20:07:19.5 


-48:22:13 


-20.95 


0.95 


G 


0.123 


Non-spherical isophotes; Intra-group emission 


NGC 6868 


20:09:54.1 


-48:22:46 


-21.21 


0.97 


G 


0.096 




NGC 6958 


20:48:42.6 


-37:59:51 


-20.68 


0.86 


C 


0.122 


Multiple shells 


NGC 7029 


21:11:52.0 


-49:17:01 


-20.26 


0.86 


G 


0.085 


Boxy inner isophotes 


NGC 7144 


21:52:42.4 


-48:15:14 


-20.51 


0.91 


G 


0.100 




NGC 7196 


22:05:54.8 


-50:07:10 


-20.47 


0.91 


G 


0.171 


Shell 


NGC 7192 


22:06:50.1 


-64:18:58 


-20.70 


0.92 


F 


0.096 


Shell 


IC 1459 


22:57:10.6 


-36:27:44 


-20.73 


0.96 


G 


0.137 


Multiple shells 


NGC 7507 


23:12:07.6 


-28:32:23 


-20.36 


0.94 


F 


0.084 


Faint N shell 



^J2000.0 

^Taken from Tully (1988) and corrected to Ho=70 km s-^ Mpc-^ 
^Colors within effective radius from Michard (2005) and de Vau- 
couleurs et al. (1991) 

Environment: C=cluster; G=group; F=field; U=unknown 
^Corrected tidal parameter as derived in subsection 3.2 
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Fig. 5. — A comparison between the initial and the cor- 
rected tidal parameter values of the sample galaxies. The vertical 
dashed and dot-dashed lines show the tidal parameter threshold for 
marginally and clearly interacting systems, respectively. The hor- 
izontal dashed line is shown for comparison with vD05 and marks 
their interaction threshold. A visual inspection of the data sepa- 
rates galaxies with detected tidal features (red circles) from visually 
appearing relaxed systems (blue triangles). 

3.2. Tidal Parameter Derivations 

Quantitative analysis of the tidal disturbance of the 
galaxies was carried out following the method described 
by Colbert et al. (2001) and vD05. Each galaxy was 
fitted with an elliptical galaxy model using the ellipse 
task in IRAF in three iterations, allowing the center po- 
sition, position angle and ellipticity to vary with radius. 
Foreground stars and background galaxies were masked 
using SExtractor^ a source detection software (Bertin & 
Arnouts 1996). Masked object frames were then divided 
by the galaxy models and were further binned by a 3x3 
factor in order to increase the signal-to-noise ratio. The 
resulting tidal parameter is given by: 



^ galaxy 



Mx,y 



(1) 



where Ix^y and M^^y are the pixel values at x,y of the 
object and model frames, respectively. 

In addition to the steps described in detail by vD05 
we applied a correction for residual noise to the derived 
tidal parameters. In order to do so we created a blank 
sky frame by aggressively masking a fully reduced object 
image and replacing the masked regions with sections 
from other object frames of the same pixel area. The 
resultant blank sky frame conserves the noise character- 
istics of the object images, consisting of both residual 
large-scale flattening variations and pixel scale photon 
noise. Each galaxy model was then added to the blank 
sky frame and a tidal parameter was re-derived for the 
model galaxy. The corrected value is then given by the 
square-root of the difference between the squared tidal- 
parameter values: 



rn2 
galaxy 



model 



(2) 



Fig. 6. — The distribution of derived tidal parameter values. The 
shaded area represents galaxies with a Tc value greater than the 
detection threshold. The four subset images are typical examples 
of various interaction levels. 

By using this correction method we allow the model 
to grow to arbitrarily large radii without erroneously in- 
creasing the tidal parameter value. This is especially 
useful when a tidal feature exists relatively far from the 
center of its host galaxy but is clearly related to it (e.g. 
NGC 596). Figure 5 demonstrates the effectiveness of 
this method which noticeably improves tidal feature pa- 
rameterization in some galaxies. 

The quantitative tidal parameter is a useful tool for 
studying the evolution of interacting galaxies. Not only 
does it correlate well with visually identified tidal dis- 
turbances, it also reveals morphological deviations from 
a smooth stellar body in visually less clear cases. Such 
an objective classification is important for comparison 
with similar surveys as well as for studying how stellar 
disturbances correlate with other galactic properties. 

3.3. Frequency of Tidal Features 

The sample of 55 galaxies can be divided into three 
categories based on their tidal parameter values. In the 
first group we include galaxies with a corrected tidal pa- 
rameter larger than 0.09. These systems all show clear 
signs of morphological disturbances both in the total- 
counts image and in the model divided frame and they 
account for roughly 53% of the sample (29 galaxies). 
The second category consists of 11 galaxies (~20%) with 
0.07 < Tc < 0.09 and it includes systems that have 
marginally detected disturbances in their stellar mor- 
phologies. Galaxies in this group typically show no ob- 
vious signs of past interactions in the total-counts im- 
age but exhibit some irregularities in the model divided 
frame. The last sub-group contains 15 galaxies (~27%) 
that lack clear interaction signatures in either the total- 
counts image or the model divided frame. All galaxies in 
this group have a tidal parameter value of less than 0.07. 

Taking Tc = 0.07 as the threshold, 40 galaxies show 
morphological signs of past interactions in their stellar 
bodies, accounting for ~73% of the sample (figure 6). 
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Fig. 7. — The distribution of tidal parameter values in respect 
to galactic environments. The red histogram represents cluster 
galaxies and the vertical lines show the mean Tc values of the dif- 
ferent populations. Cluster galaxies (Tc = 0.056) are on average 
less tidally disturbed than group and field galaxies (Tc = 0.102). 

This result confirms the findings by vD05 and further 
generalizes the phenomenon due to the statistical com- 
pleteness of this study. 

An additional sub-population of the sample consists of 
galaxies that have a projected close companion. We iden- 
tify 7 such pairs, all showing morphological disturbances 
in both galaxies. A visual inspection of the companions 
finds that 5 of them have spiral features and are likely 
disk-systems while the other 3 are probably early-type 
or dwarf galaxies with no significant disk component. In 
addition to these interacting systems we identify one pair 
(NGC 5638) that shows no interaction signatures in ei- 
ther of the galaxies and is appropriately assigned a tidal 
parameter value of 0.04. The pair fraction is similar to 
that found by vD05 and although the fraction of early- 
type companions is lower (see Whitaker & van Dokkum 
2008), the difference is insignificant. 

4. PROPERTIES OF TIDALLY DISTURBED GALAXIES 

The derived tidal parameter allows us to quantify the 
disturbance caused to the stellar body of galaxies due 
to gravitational interactions. One of the key advantages 
of such a numerical parameter is that it allows us to 
correlate gravitational disturbance signatures with other 
galactic properties. In this section we compare the de- 
rived Tc of each galaxy with galactic environment, its 
broadband optical colors and its nuclear activity. These 
properties have all been hypothesized to be related to 
gravitational interactions and the subsequent evolution 
of these galaxies. 

4.1. Galaxy Environment 

The fraction of galaxies with Tc > 0.07 in our sample 
increases to ~80% when cluster members are excluded. 



This is depicted in figure 7, which shows the relation 
between the environment in which the galaxies reside and 
their tidal parameter values. A one-sided Mann- Whitney 
test finds a 1.4% chance that cluster galaxies have a tidal 
parameter value larger than the sample median. 

Interaction signatures seem therefore to be typical of 
group and field galaxies, and relatively rare in clusters. 
The absence of tidal features in clusters might have been 
expected given that mergers are rare in relaxed clusters 
(Makino & Hut 1997; van Dokkum et al. 1999); on the 
other hand, "galaxy harassment" due to frequent long- 
range interactions has been shown to create long, faint 
tidal tails associated with infalling galaxies (e.g. Moore 
et al. 1996, 1998). 

4.2. Broadband Optical Colors 

The relation between the optical colors of elliptical 
galaxies and their absolute magnitude has been known 
for more than three decades (e.g. Sandage & Visvanathan 
1978; Frogel et al. 1978). This relation states that giant 
ellipticals have systematically redder colors than their 
less massive counterparts. In recent years this observed 
phenomenon has been mostly attributed to the evolution 
of these galaxies through hierarchical growth and mass 
accretion, thus differentiating it from monolithic collapse 
models. This change of paradigm has been supported by 
both observations and numerical simulations of galactic 
evolution (e.g. Kauffmann & Chariot 1998; van Dokkum 
et al. 1998; Shioya & Bekki 1998). 

In their paper, Schweizer & Seitzer (1992) show that 
the observed "fine structure" value of ellipticals is also 
correlated with the systemic colors of these galaxies. The 
correlation persists even after the observed values are 
normalized by the known color-magnitude relation. We 
perform a similar analysis of our data using normalized 
optical colors in conjunction with the objectively derived 
tidal parameter (see subsection 3.2 for details). 

Figure 8 shows the color-magnitude relation of the 55 
sample galaxies for the B-V color, plotted against ab- 
solute V magnitudes. Color measurements were taken 
from Michard (2005) where possible (41 galaxies) and 
from de Vaucouleurs et al. (1991). This relation has an 
inherit la scatter of 0.042 in color which cannot be all at- 
tributed to errors in measurement and it agrees with the 
slope found by Schweizer & Seitzer (1992). We subtract 
the color-magnitude relation from the measured colors 
and plot the residuals against the corrected tidal param- 
eter values of the galaxies. Figure 9 shows that tidally 
disturbed galaxies are on average bluer than their non- 
interacting counterparts, in agreement with the findings 
of Schweizer & Seitzer (1992). This correlation is moder- 
ate, showing a difference of only 0.05-0.1 in color between 
galaxies with tidal parameter values of 0.0 and 0.2. A 
Spearman rank test finds that the probability that this 
correlation was drawn by chance is less than 0.01%. 

The correlation between tidal disturbances and broad- 
band colors suggests that these merger events probably 
trigger some star formation in the interacting galaxies. 
Alternatively, it is possible that galaxies accreted into el- 
lipticals typically contain a significant fraction of young 
or metal-poor stars. Such a scenario is not uncommon in 
groups, where a single massive elliptical is surrounded by 
multiple low mass late- type galaxies. Nevertheless, it is 
important to note that all the galaxies in the sample are 
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Fig. 8. — B — V color-magnitude relation of the sample galaxies. 
The solid line is a fit to the data and the dotted lines represent la 
deviations from the fit. 

red and that any contribution to their broadband colors 
from star- formation or young steUar populations is small. 

4.3. Tidal Features and the AGN Duty Cycle 

It has been suggested by several authors that gravita- 
tional interactions play an important role in triggering 
galactic nuclear activity. Some of these studies show a 
correlation between the brightness of elliptical galaxies 
in radio continuum observations and morphological dis- 
turbances (e.g. Heckman et al. 1986; Smith & Heckman 
1989). In this scenario minor merger events induce dust 
and gas accretion onto the galactic center and then trig- 
ger nuclear activity. The same mergers should in princi- 
pal leave their mark on the stellar body of the accreting 
elliptical in the form of morphological disturbances such 
as described in subsection 3.1. 

In order to test this hypothesis we compiled a list of 
6cm and 20cm flux measurements for 26 galaxies with 
publicly available data. The data were collected form 
a number of catalogs, including the NRAO VLA Sky 
Survey (Condon et al. 1998), the FIRST survey (White 
et al. 1997), the N0RTH6CM and the NORTH20CM 
databases (Becker et al. 1991; White & Becker 1992), 
the Green-Bank 6cm survey (Gregory et al. 1996) and 
the Dixon Master List of Radio Sources (Dixon 1970). 
As not all galaxies have published measurements for both 
wavelengths we derive a correction factor from the 20cm 
to 6cm relation using 8 galaxies for which both data are 
available. This relation is expressed by the following 
equation: 



logF6cm = 1.20 + 0.53 log F20cm 



(3) 



where Fq^u and F20m are the measured fluxes in the 6cm 
and 20cm bands, respectively. The la scatter in this 
relation is 0.29 dex. Figure 10 shows the correlation be- 
tween radio continuum flux and tidal parameter for the 
sample galaxies. Plotted values are either 6cm observa- 
tions or 20cm observations corrected to 6cm using the 
linear relation from equation 3. 
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Fig. 9. — The B-V color is normalized by the color-magnitude 
relation from figure 8 and plotted against the corrected tidal pa- 
rameter. There exists a moderate correlation between tidal dis- 
turbances and the broadband colors of luminous elliptical galaxies. 
The solid line is a fit to the data and the dotted lines represent la 
deviations from the fit. 

We see no relation between the two parameters. This 
lack of a tight correlation is itself not surprising given 
the very different expected lifetimes of tidal features and 
the radio duty cycle. Tidal interactions, as portrayed 
by their assigned tidal parameter values, are expected 
to survive for a long time in the outskirts of elliptical 
galaxies. This means that a disturbed galaxy can retain 
its large value of tidal parameter even after the nuclear 
activity has quieted. 

This can be a viable explanation for the scatter at high 
values of Tc but not for that at the low end of the Tc dis- 
tribution. Since the lifetime of the radio mode is only 
roughly 10^ years (e.g., Croton et al. 2006; Shabala et al. 
2008), galaxies with a low Tq value are all expected to 
be quiet at radio wavelength. The existence of radio- 
loud AGN in undisturbed ellipticals is very interesting 
and may imply that gravitational interactions are not 
the only and possibly not the most important AGN trig- 
gering mechanism. 

5. ACCRETION RATE AND INFERRED MASS GROWTH 

In the following section we will attempt to estimate the 
rate at which ellipticals grow through mergers using two 
techniques. First we estimate the survival time of the 
observed tidal features using the dynamical time of the 
system at the innermost radius of feature detection. This 
estimate relies on the assumption that tidal features dis- 
sipate due to dynamical mixing in their environment. As 
an alternative approach we analyze the interacting pairs 
in the sample using friction timescale approximations as 
suggested by, e.g., Patton et al. (2000). 

5.1. Typical Tidal Feature Survival Time 

Most of the merger rate estimates in the literature 
measure the frequency and properties of close interact- 
ing pairs. In this work, however, we are able to use an 
alternative approach. For each galaxy we measure the 
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Fig. 10. — The relation between radio continuum flux and tidal 
parameter for 26 of the sample galaxies. The existence of radio- 
load, non-interacting galaxies suggests that gravitational interac- 
tions are not the only, and perhaps not the most important, AGN 
triggering mechanism. 

smallest radius where tidal features are seen and esti- 
mate the age of the merger event from the dynamical 
time at that radius. 

The dynamical time at a distance r from the center of 
a galaxy with circular velocity Vc{r) is defined as the time 
it would take a test particle to complete one-quarter of an 
orbit under the gravitational potential of a homogeneous 
sphere. We use this to assess the survival time of a tidal 
feature, assuming that after completing No orbits the 
feature will dissolve into the smooth stellar distribution 
due to the inherent velocity dispersion and differential 
rotation of the galactic potential. The survival time is 
then given by: 



Ts = NoTdyn — No 



2 '^c('^min) 



(4) 



where No is the number of dynamical times it takes tidal 
features to dissolve into the stellar body of the galaxy and 
Vmin is the smallest radius at which we detect tidal fea- 
tures. In order to find this radius we derived fractional 
tidal parameter values as a function of radius for the 
sample galaxies, thus measuring tidal distortions within 
limiting radii. Figure 11 shows the average radius of each 
fractional Tc bin of size ATc = 0.005. The majority of 
the significant tidal features lie at radii larger than 23 
kpc which is also the average distance of features with 
Tc ~ 0.07. We estimate Vc by measuring the circular 
velocity of a model galaxy with an NFW dark matter 
profile (Navarro et al. 1997) and a Hernquist stellar pro- 
file (Hernquist 1990) at the same radius. The typical 
time since the last accretion event for these galaxies is 
then: 



Ts = 0.35No 



23 kpc 



395 km s 



Gyr, (5) 



The actual survival time, however, can vary greatly 
based on the tidal feature morphology and kinemat- 
ics that result from different interactions and collisions. 



Since typical stellar orbits in elliptical galaxies are highly 
eccentric and have significant precession, radial features 
will dissolve much quicker than tangential ones. For 
example, a radial linear feature is only expected to be 
detectable for about one dynamical time {No ~ 1), 
whereas shell-like features are expected to survive for 
many {No > 10). 

This technique of estimating the interaction time-scale 
of nearby ellipticals is important as it does not rely on 
close pairs, which are much less common. Nevertheless, 
when deriving the survival time we did not treat any 
biases that might be caused by the projection of tidal 
features onto the plane of the sky. For example, a stellar 
ring around a galaxy at some distance r may appear as a 
radial feature extending all the way through the center of 
the galaxy if aligned properly. Rather than individually 
simulating all sample galaxies we assume that feature ori- 
entations are random and that our sample is sufficiently 
large to minimize their effect on this result. 

5.2. Implied Merger and Mass Growth Rate 

Assuming that in a given merger event all the mass 
from both galaxies is retained in the system we can write 
the mass growth rate of an accreting elliptical: 



dM _ m2 fr 
M ~ miTs 



(6) 



where m2/mi is the progenitor galaxy mass ratio, Ts is 
the tidal feature survival time (see equation 5) and /t is 
the fraction of observed galaxies with tidal distortions. In 
this paper we confirm the value of /t ~ 0.7 and measure 

-J- s\' min 

) - 0.35Aro Gyrs. 
Two significant uncertainties remain in equation 6. 
Both the progenitor mass ratio m2 / mi and the number of 
survival orbits No cannot be observed directly from our 
data and require simulations of multiple systems in order 
to constrain their values. In addition, both parameters 
can likely have a wide range of values, varying between 
systems and tidal feature morphologies. We therefore 
adopt conservative values for these parameters to get an 
estimate for the mass growth rate due to gravitational 
interactions: 



dM 



0.2 



1712 /mi 
0.1 



A 

0.7 



dT 



(7) 



Assuming a constant rate since z ~ 2, this implies 
that an isolated galaxy can grow up to 3 times its initial 
collapsed mass by minor mergers only. This result is in 
good agreement with recent studies that show that galax- 
ies can grow significantly in mass through minor mergers 
(Naab et al. 2007; Bournaud et al. 2007; Bezanson et al. 
2009). 

5.3. Pair Friction Time 

For purposes of comparison with previous studies we 
perform a second analysis of the interaction timescale 
using a subsample of close pairs of interacting galaxies. 
In these systems the time it would take the two galaxies 
to merge can by approximated by the dynamical friction 
timescale (Patton et al. 2000): 



2.64 X lO^r^Vc ^ 
= MlnA 



(8) 
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Fig. 11. — The average distance at which tidal features are 
detected. The derived fractional tidal parameter values are binned 
with bin size of 0.005. The horizontal line marks our detection 
threshold of Tc=0.07. Most detections are to the right of the 23 
kpc line which is the average at the 0.07 threshold. 

where r is the physical separation of the pairs, Vc is 
the circular velocity, M is the mass of the secondary 
galaxy and A is the usual Coulomb logarithm (Binney 
& Tremaine 1987). In our galaxy sample we identify 7 
merging pairs with a mean projected distance of 21.4zb6.7 
kpc, resulting in a circular velocity of '^400km s"""^. Us- 
ing typical values of Mc = 7 x IO^^Mq and InA = 2 
we get a second interaction time estimate of ^0.53 Gyr. 
This result is in good agreement with vD05 and it is 
similar to Tg for Nq ^ 2 

6. CONCLUSIONS 

The evolution and growth of elliptical galaxies in the 
nearby Universe is driven by gravitational interactions 
and merger events. The rate of these interactions, which 
are not accompanied by significant star formation, can be 
estimated from the morphological distortions that they 
cause to the stellar bodies of colliding galaxies. vD05 
showed that at z^O.l signatures of past interactions are 
observed in most ellipticals and can be accounted for in 
more than 73% of the galaxies. In this work we con- 
firmed this result upon analyzing 55 nearby ellipticals in 
a variety of galactic and environmental properties. This 
is the first time that such a study is performed on a sta- 
tistically complete sample and it shows, above all, that 



elliptical galaxies are still assembling via mergers and 
accretion events in the nearby Universe. Although this 
high rate of gravitational interaction signatures is not a 
new finding, the completeness of our sample implies that 
this is an important phenomenon for all ellipticals and 
that it is not restricted to certain populations. Moreover, 
the high frequency of tidal features strongly suggests that 
these interactions play a critical role in the evolution and 
growth of elliptical galaxies. 

In this paper we explored the link between tidal dis- 
turbances and other evolution driven properties of the 
sample galaxies. We found that interacting systems 
have slightly bluer broadband colors than non interacting 
ones, implying that these mergers are accompanied by lit- 
tle or no star formation. One explanation for this is that 
the progenitor galaxies are mainly composed of old stellar 
populations, thus producing a red galaxy. Alternatively, 
the total flux from young stars in the accreted galaxies is 
small compared to the accumulated light from old stars 
in the ellipticals, suggesting that these are mainly minor 
merger events. 

We also examined the proposed link between gravita- 
tional interactions and nuclear activity in ellipticals (sub- 
section 4.3). We accomplished this by relating the radio 
continuum flux of the sample galaxies to their respective 
tidal parameter values. Given the long dynamical time 
in the outskirts of ellipticals, tidal features, unlike radio 
activity, are expected to survive there for a long time. 
Using equation 5 it is easy to show that at a radius of 
50kpc a tidal feature will remain for ^ 10^ years. On the 
other hand, the standard assumptions predict a lifetime 
of only 10^ years for the AGN radio mode. This suggests 
that all radio loud galaxies are expected to have a high 
value of tidal parameter, in contradiction to our findings 
(figure 10). We conclude that gravitational interactions 
are therefore an unlikely candidate for provoking nuclear 
activity in elliptical galaxies. 

Lastly, we estimated the lifetime of observed tidal fea- 
tures by measuring the dynamical time of the galaxies at 
the innermost radius of feature detection. We found that 
the sample galaxies are generally relaxed at radii smaller 
than ^23 kpc, implying that nearby ellipticals typically 
experience gravitational interactions every ~0.35 Gyrs. 
This value yields a mass accretion rate of dM/M~0.2 
with large uncertainty. The derived rate of mass growth 
therefore shows that elliptical galaxies can grow signifi- 
cantly through (presumably mostly minor) mergers and 
low-mass accretion events. 



APPENDIX 

In this section we present the full data atlas. All the images were smoothed using a box-car kernel of dimensions 
8" x8" and were further bin by a factor 2x2 in order to bring out some of the fainter tidal features. Similar to figure 
4, every line shows from left to right the full data set, object masked frame, the fitted model and the model-divided 
residual image. The corrected tidal parameter and a scale bar are also shown. 

The full paper and appendix are available online at: http://www.astro.yale.edu/obey 
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